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The question of how fast a quantum state can evolve has attracted 
a considerable attention in connection with quantum measvue- 
ment, metrology, and information processing. Since only orthog- 
onal states can be unambiguously distinguished, a transition from 
a state to an orthogonal one can be taken as the elementary step 
of a computational process/ Therefore, such a transition can be 
interpreted as the operation of "flipping a qubit", and the num- 
ber of orthogonal states visited by the system per unit time can 
be viewed as the maximum rate of operation. 

A lower bound on the orthogonalization time, based on the 
energy spread AS, was found by Mandelstam and Tamm.^ An- 
other bound, based on the average energy E, was established 
by Margolus and Levitin.^ The bounds coincide, and can be ex- 
actly attained by certain initial states if AE = E. However, the 
problem remained open of what the situation is when AE 7^ E. 

Here we consider the unified bound that takes into account 
both AE and E. We prove that there exist no initial states that 
saturate the bound if AE 7^ E. However, the bound remains 
tight: for any given values of AE and E, there exists a one- 
parameter family of initial states that can approach the bound 
arbitrarily close when the parameter approaches its limit value. 
The relation between the largest energy level, the average energy, 
and the orthogonalization time is also discussed. These results 
establish the fundamental quantum limit on the rate of operation 
of any information-processing system. 

Starting with the classical result of Mandelstam and Tamm,^ 
it was later shown by Fleming,* Anandan and Aharonov,^ and 
Vaidman® that the minimum time r required for arriving to an 
orthogonal state is bounded by 



r > h/iAE, 



(1) 



where {AEf = WH^W - {{ip\H\ip)f, H is the Hamiltonian, 
and the wavefunction of the system. A different bound was 
obtained in,^ namely, 

T > h/AE. (2) 

Here, E = {%l)\H\'ip) is the quantum-mechanical average energy 
of the system (the energy of the ground state is taken to be 
zero). Both bounds (O and ^ are tight, and achieved for a 
quantum state such that AE = E. 

Since then, a vast literature has been devoted to various as- 
pects of this problem. In particular, inequality ([2} has been 
proved for mixed states and for composite systems both in sep- 
arable and in entangled states (e.g., Giovannetti et al.,'^'* Zan- 
der et al.^). Bound ([2]) obtained for an isolated system has 
been generalized to a system driven by an external Hamilto- 
nian (a "quantum gate") in.^"''^^ Various derivations of ^ 
and ([2]) (e.g.,^^"^**), bounds based on energy-distribution mo- 
ments, more general problems of time-optimal quantum evo- 
lution, and the ultimate limits of computation'^^' '^^ have 
been considered. 

However, what remained unnoticed is the paradoxical situa- 
tion of the existence of two bounds based on two different char- 
acteristics of the quantum state, seemingly independent of one 
another. Since the average energy E and the energy uncertainty 
AE play the most determinative role in quantum evolution, it is 
important to have a unified bound that would take into account 
both of these characteristics. 



In all known cases where bounds ([1} and ((2)) can be exactly 
attained, the ratio a = equals 1. A question arises: what 
happens if a 7^ 1? Some authors just assumed without justifi- 
cation that the minimum orthogonalization time is 

r„>i„ = max ( — , ^) = ^{E + AE ^ \E - AE\- 

In fact, the situation is not so simple. Bound indeed, can 
only be achieved for a = 1. However, this bound remains tight 
when Q 7^ 1 as well, though in this case it is only asymptotically 
attainable. 

Theorem 1 Under the assumption that the smallest (ground) 
energy of a quantum system is zero, 

1. The only state that attains bound {!]) is the two-level state 

W^i^Wo) + \M), (4) 

where H\ijjh) ~ kEi for k = 0, 1; 

2. The only state that attains bound {2) is likewise state {3)). 
State {4)) is unique up to degeneracy of the energy level E\ and 
arbitrary phase factors for It/jq) and IV'i)- 

Proof. To prove the first statement we shall use the 
trigonometric inequality 



cos 2: > 1 ^ a:: sin a; -x^ , (5) 



t2' 



which is valid for any real x. Note that turns into an equality 
iff a:: = or a; = ±7r. 
Let the initial state be 



|V'(0)> = ^c„|i5„), 



(6) 



where the \En) are energy eigenstates of the system and 
Er=o|cnP = l- Then 



= E 

n.n' 
00 
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(7) 
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Using inequality ([5}, we obtain 



h/t 



\S{t)\'>l~- J2 1=' 
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h/t 



(8) 



U d\S{t)\'^ 1 / AEy 

Since \S{t)\'^ > 0, it follows that ^^^^ = whenever S{t) = 
0. Thus, at a time r such that 5(r) = 0, the second term in 
(151) vanishes, and we obtain 



> 1 



ill 

7r2ft2 



(AEf 



(9) 



1 



which yields inequality ([T}; this is just another way to derive 
that bound. However, for ^ to turn into an equality, it is 
necessary that inequality ^ turn into an equality for every 
term of the double summation ((Tjl . Hence, 



either x„„i = 



En — E^i 
h/T 



En — E„i 



0^ n -^n' I 
or x„„f = — = ±7r 

h/T 

(10) 

for all n, n' such that Cn 7^ 0, c„' 7^ 0. It follows that, to attain 
bound IJ]), |^/'(0)) must be a superposition of only two energy 
eigenstates with energies Eo = Q and Ei. 

To prove the second statement, we repeat briefly the deriva- 
tion given in.^ This time we use the trigonometric inequality 



cosx > 1 (a; + sinx). 



(11) 



valid for all x > 0. Again, turns into an inequality only for 
a; = or X = TT. 
Then 



Ent 



ReS(t) = |c„|^cos 

>f:i-i^[i-f(^+sin 



h 

2 / Ent 



Ent- 



(12) 



n = 

^l^^ + hmS{t). 

At time r. Re S(r) = Im S(r) = 0. Hence > 1 - 2|* , which 
results in bound ([2]). However, to actually attain this bound, 
inequality (lllf) must become an equality for every term of the 
sum (|12p : that is, for any n, 

Ent Ent 

either Xn ~ — ^ = or Xn = — ~ vr; 

n ft 

This is possible iff \i^{0)} has the form (jU. 

Let us show now that no mixed state can attain bound ((3)1. 
Consider an initial mixed state density matrix p(0) with spec- 
tral decomposition 



p(0)=^A,|V«(0)>(^«(0)|, 



(13) 



where Ai > 0, J2i>^i = 1. and (V'^'Ho)!^*^' (0)) = At time 
t, the density matrix becomes 

p(t) = ^A.|^«(t))(V'<')(t)|. 

i 

Hence, if r is the orthogonalization time, then 

Tr[p(0)p(r)] = ^^A,A,l(^W(0)|^(^-)(t))|^ = 0. (14) 



Note that all terms in the above sum are non-negative, and 
therefore all of them must be zero to satisfy (|14p . As shown 
above, terms (V'''>(0)|V''*'(t)) = 0, for T given by ®, iff each 
Itp^'Ho)) has the form ©. 

Consider two functions from spectral decomposition (|13p 



^'^'(0)) = 75(1^0) + |^r'(0))) and 



V2 



(15) 



where (t/''''(0)|V'*^'(0)) = and |a|^ + \bf = 1 (the zero-energy 
state is nondegenerate) . Since 

(^«(0)|^(^)(0)> = i((V.o|v^o) + a(^«(0)|Vf^(0)» =0, 



it follows that a = -1 and \ip^'\0)} = -^{\i^{0)) ~ \->pi"Ho))). 
But then 

{^W(0)|^(^'(r)) = i((V>o|^o) - (^f'(0)i^«(r))) = i / 0. 

Hence, equality (|14|l cannot be satisfied for r given by 
Thus, this bound is not attainable by a mixed state. ■ 

For pure states, results of Theorem[T] also follow from the 
analyses presented in^'^ and.^^ 

Theorem[l] shows that bounds ([1} and ((2)) can only be at- 
tained by a state for which a = = 1. Such a state is unique 
up to degeneracy of the energy level Ei and arbitrary phase 
factors for l^po) and l-i/Ji). It follows from that theorem that 
there exists no initial state \ip{0)) with a 7^ 1 that would attain 
bound ((3]), i.e., no state can evolve into an orthogonal state in 
the minimum time given by ([T]) or ((5)). The question to be an- 
swered is how close it is possible to approach the unified bound 
([3]). Let us rewrite (|3]) in a different form. 



/ h 

'mm — max , 

V 4:E 



h \ _ fe(l + el'"°l) 
IAe) ^ AEil + a) 



(16) 



Theorem 2 For any a ^ 1 and any e > 0, there exists a state 
\-ipiO)) such that (V'(O)IV'(r)) = at time 



r < 



4E{l+a) 



'-(l+e). 



(17) 



Proof. We will show that in both cases, a < 1 and a > 1, 
there exist families of initial states that approach limit (|16|) 
arbitrarily close. 

1. Let a < 1. Consider a family of states 

[^^(O)) = co\0) +ci\Ei} +C2\E2). (18) 
Denote |cip ~ Pi, and introduce dimensionless variables 

x,{t) ^ 2tvE,^, 0,1,2 (19) 

(the rotation angles of the state vectors \Ei)). Then S{t) — 
{iP{0)\^{t)) = iff Xi{r) = Xi such that 



and 

Of course. 



pi sin xi_ + p2 sin X2 = 

PO + Pi cos Xi -I- P2 cos X2 = 0. 

Po +Pi +P2 = 1. 



(20) 

(21) 
(22) 



Let < po << 1 be a parameter of the family of states (|18|l 
(one cannot set po = 0, since that would change the ground 
energy level). When po is small, the values of xi and X2 differ 
almost exactly by tt. Let 

X2 = TT -|- xi — 5sinxi, with (5 << 1. (23) 

Then, from ((20)) . 

Pi =P2[1 -(Jcosxi -|-0((5^)]. (24) 

Substituting into (PT)) and (f^ we obtain 

pi = i - |(H- cos xi) -1-0(5''), 

P2 = |-|(l-cosxi)-FO(<5'), (25) 

P0 = l + O{5''). 



Using (|23p and (|25p . we can calculate the average value {x) — Thus, choosing a sufRciently large k > guarantees that 



-E, namely, 

(x) = piXi +P2X2 

= ^+xi~ |[2a;i - 2sina;i +7r(l - cosa;i)] + 0(5^). 



The standard deviation Ax = ^^AE is obtained as 



Ax = ^Jpixl +P2xl - (piXx +P2X2Y 

= f + ^[a::i(7r + xi) - vr sin xi] + 0(5^ 
On the other hand. 

Ax TT 



Hence, 



E 



Xl 



(x) n + 2xi 



0(5). 



(26) 



(27) 



(28) 



0(5). 



Substituting (|29l) into (|27[l and taking into account the last 
expression from (|25|l yields 

Aa: = I + po [f - 1) - sin f (i - l)] + 0(pg), (30) 

which results in 
hAx 



2ttAE 
4AE L 2 Va2 



- 1 - 



lsin2^(i-l))+0(pg)]. 



(31) 



Finally, choosing 

po < 2e [-!^ - 1 - - sin f (i - l) 
we obtain 

h ^ h . 

4AE<^^4AE^' + '^- 
2. Let a > 1. Consider a family of states 

|V(0)) = Co|0) +Ci|£i) +C2fc+l|S2fc+l), 

where 

fe = 1, 2, . . . ; E2k+i = (2fe + l)Ei- \cof = po = 

|C1 1^ = Pl = 1(1 - ■^) ; |c2fc+l 1^ = P2fc+1 = jfj- 



(32) 
(33) 

(34) 



Using the dimensionless variables p9[) . it is readily seen from 
(|33p and (|34|) that S{t) = at the least time r for which 



h ^ h , 

< r < (1 + e). 

AE - AE^ ' 



(40) 



This completes the proof. ■ 

As graphically illustrated in Fig.[T] Theorem[5] means that 
there are families of states that approach equality in ((3)l for a 
given value of a, in the limit as a certain parameter approaches 
zero (for a < 1) or infinity (for q > 1), but no state yields a 
strict equality in ([Sjl for a 7^ 1. It turns out that bound ([T]) is 
tight whenever a = < 1, while ([2]) is tight whenever a > 1. 



(29) 3 








3 



Figure 1: The solid line shows bound (|16|l : the dotted lines — 
parametrized, respectively, by po — > and k — > 00 — show suc- 
cessive approximations to this bound, corresponding to initial 
states from families (1181) and (I33II. 



Another interesting question is the relationship between the 
maximum energy eigenvalue -Emax that contributes to \^I){Q)), 
the orthogonalization time, r, and the average energy E. Note 
that determining the minimum r for a given E is equivalent to 
determining the minimum E for a given t. The next theorem 
provides simple but useful results. 

Theorem 3 Let r be the minimum time such that S{t) — 
{tp{0)\tp{T)) = 0. Then there exist a state |'i/'(0)) and a set 
of energy eigenvalues {E„} of all energy eigenfunctions that 
contribute to \tpo) such that 



E„ 



h 



and 



E„ 



xi{t) = Xl = ^E^ = TT and X2k+i = ^^{2k+l)E^ = 7r(2fc+l). 
Then 



<E< 



(x) = piXi + p2fc+ia;2fe+i = f (1 + x) 



(35) 
(36) 



Proof. Let 



and 



S(r) = ^|c„|^e-^^'^"^/''=0. 



(41) 
(42) 

(43) 



Ax = \JPlxl +P2k+ix\^+i - {piXl +P2k+lX2k + l) 

= f(l + 4/3+f) + 0(^). 
Hence, a = = 1 + 4/3 + 0{\) and 



Since (x) = '^E, it follows from (|3Z| and JSlJ that 



(37) 



(38) 



(39) 



Suppose Ek > 7. Then = £fc - | > and e'^'^'-^''^/'' = 
g-2^i£;^^r/h Ylence, replacing Ek by E^ in (|43| will not affect 
the equality. Thus, the same ortogonalization time r can be 
achieved with smaller average energy E' = E~ \ck\^h/T. This 
proves inequality (|4ip . 

Now, let -Bmax be the largest energy of the energy eigenfunc- 
tions that contribute to |'i/>(0)), and let the average energy be 
iJ^^'. Obviously, the validity of 1)43^ will not be affected if we 
replace all energy levels _E„ by Smax — En- Then the average 
energy will become S'^' = E 

max Since E can be chosen 



as min(£;(i',£;(2'), this proves that E < Smax/S. Also, substi- 
tuting (HI) into (|IT1) yields i?max/4 < E. Thus, (HH) is proved. 
■ 

Theoreni[3] allows us to restrict our attention to states that 
satisfy (|4l]l and (|42|) for further analysis of the orthogonal- 
ization time r. 

Theorem 4 Let E'max be the maximum energy of all en- 
ergy eigenfunctions that contribute to \tp{0)). Then the 
minimum orthogonalization time t obeys the inequality 



T > 



2E„ 



(44) 



and equality is attained iff |'0(O)} has the form HP with 



Proof. Let 



S{r) = immr)) = J2 \cn 



\1 -l-KiE^lh _ 



0, (45) 



n=0 



where i?„ 



E„ 



Then, if T < 



2E^ 



all terms in the 



above sum except the zeroth have strictly positive imag- 
inary parts — which violates equality (|45p . Similarly, if 
r = 2e!^' then jS]) is satisfied iff |^(0)) is of form (g]) 

with El = L^max- ■ 

The equality case in (|44| had been considered in.^'' 

Some authors (e.g.,^"*) prefer to write expressions ([2]) and 
(glD as r > h/^{E-Eo) and r > h/2{E^^^- Eq), arguing 
that a constant shift of the energy spectrum affects only 
the overall time-dependent phase of the wavcfunction. In 
our opinion, however, it is just this "freedom of shift" that 
allows one to set to zero the energy of the ground state 
even if the smallest eigenvalue of the Hamiltonian is not 
zero (just as in the case of the hydrogen atom). 

The authors wish to thank Prof. D. Brody (Imperial 
College, London) for a useful discussion. 
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